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ABSTRACT 

It is thought that planetary mass companions may form through gravitational 
disk instabilities or core accretion. Identifying such objects in the process of 
formation would provide the most direct test for the competing formation the¬ 
ories. One of the most promising candidates for a planetary mass object still 
in formation is the third object in the FWTau system. We here present ALMA 
cycle 1 observations conhrming the recently published 1.3 mm detection of a dust 
disk around this third object and present for the first time a clear detection of a 
single peak (2-1) line, providing direct evidence for the simultaneous exis¬ 
tence of a gas disk. We perform radiative transfer modeling of the third object 
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in FW Tan and find that current observations are consistent with either a brown 
dwarf embedded in an edge-on disk or a planet embedded in a low inclination 
disk, which is externally irradiated by the binary companion. Further observa¬ 
tions with ALMA, aiming for high SNR detections of non contaminated gas lines, 
are required to conclusively unveil the nature of the third object in FWTau. 

Subject headings: circumstellar matter — protoplanetary disks — stars: individ¬ 
ual (FW Tan) — planetary systems — techniques: interferometric 


Introduction 


Planetary mass companions (PMC) at wide separations often exceeding 100 an have 


been relatively frequently found in recent direct imaging surveys (e.g. 

Neuhauser et al.||2005 

Lafreniere et al.|2008 

Schmidt et al. 

2008 

Ireland et aL|2011 

Bailey et al.j2014 

). If relatively 


close, i.e. at separations < 100 an, PMCs may form through the classical mechanisms of core 
accretion (Lissauer & Stevenson 2007) or Class II gravitational instabilities (Boss 1997, 2011[). 


At larger separations, disk fragmentation at the Class O/I stages seems to be more likely 
(Kratter et al.||20l0). Due to their large separations from their host stars, these objects are 


probably the most promising candidates to directly observe and characterize circumplanetary 
disks that may have formed either from the surrounding cloud or from the massive disk 
around the host star. Indirect evidence for the presence of such disks around several PMCs 
has been provided recently by the detection of emission lines or possible mid-infrared excesses 


that might be related to disks and outflows (e.g. Seifahrt et al.|2007 Schmidt et ah 2008 


Bowler et ah 2014), and large optical/ultraviolet excess emission, potentially indicating the 


presence of shocks produced during the accretion process (Zhou et ah 2014). 


One of the most convincing candidates for a disk around a potential PMC has been re¬ 
cently provided by the third component in the triple system FW Tan. This system belongs to 


the nearby {d ~ 140pc) and young (~ 2Myr) Taurus-Auriga star forming region (Goldsmith 


et ah 2008). The primary in FWTau consists of a close binary, composed of two M5 stars. 


with a projected separation of 11 an (75mas). The faint tertiary component has been hrst 


reported by White & Chez (2001) and conhrmed to be co-moving at a projected separation 


of 330an by Kraus et ah (2014). The primary binary in FWTau is not showing evidence 


for the presence of accretion (Cieza et ah 2012) while for the third object indications for 


accretion have been found (Bowler et ah 2014|). 


Moreover, Kraus et al. (2015, hereafter K2015) have recently derived a dust disk mass of 


1-2 Mm, based on their ALMA band 6 continuum observations. While their observations also 
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cover the (2-1) line, their observational setup used wide velocity channels (20knis“^) 
and the line remained undetected. Despite these constraints however, the nature of the 
tertiary companion in the FW Tau system has remained unclear, as the near-IR photometry 
and spectroscopy is consistent with both a planetary mass object with accretion-induced 
veiling, or a brown dwarf or low-mass star with spectral types M5-M8 embedded in an 
edge-on disk (Bowler et ah 2014). 


Here we present a clear detection of the (2-1) line towards the tertiary companion 
to FW Tau and construct the spectral energy distribution (SED) of the system from near- 
IR to the millimeter wavelengths. Using radiative transfer and CO-line disk modeling we 
explore the two proposed scenarios. 


2. ALMA cycle 1 observations: continuum and gas detection 


We observed FWTau with ALMA in Cycle 1 on 2013 December 4. We obtained one 
epoch observation in band 6, with the correlator conhgured to obtain one baseband centered 
at 230.52 GHz which was aimed at detecting the (2-1) spectral line, and three contin¬ 
uum basebands centered at 228.52, 214.52, and 212.52 GHz. The total bandwidth for the 
observations was 7.5 GHz, with a unique spectral spacing of 488.28 kHz in 3 840 channels for 
each 1.875 GHz baseband. 


The observations were carried out using 27 antennas in a compact configuration, im¬ 
plying we obtained baselines ranging from 15.8m to 462.9m (12.2 to 356.1 kA). Standard 
calibration steps were applied to the data. The calibration sources associated to these ob¬ 
servations were QSO J042315-012033 for band-pass calibration, and QSO J051002-I-180041 
for gain and phase calibrations, and one antenna was flagged-out because of high system 
temperature. The observations consisted of 4 scans, which translates to a total time on 
source of 3.6 min for the field. After calibrating the set of raw visibilities, we applied one 
iteration of self-calibration in order to correct low-order phase calibration errors. 


Deconvolving the set of visibilities with the GLEAN task implemented in GASA (Mc- 


Mullin et al.||2007) and using natural weighting we obtained the final images. The continuum 


image has an rms of 0.18mJy/beam, with an elongated beam size of 0'.'75 x 1731 and a po¬ 
sition angle of —48.6° (north-east). The mean rms associated to the individual ^^GO (2-1) 
channels is 21.1mJy/beam per 0.635 kms“^ bin, after removing the continuum contribution 
in the visibility domain. 


The continuum image shows a clear detection centered on the location of the faint 
companion, while it shows no significant emission from the binary system, confirming the 
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detection reported in K2015 The total flux from the third object is 3.4 ± 0.2 mJy. This 


detection is above the 1.78 ± 0.03 mJy reported by K2015 and in marginal agreement with 
the 4.5 ± 1.1 mJy at 850/xm measurement of Andrews & Williams (2005). 


Similar to the continuum observations, we detected weak emission from (2-1) 

gas located at the same position as the third object, with a SNR ~ 4 in three consecutive 
channels. This emission appears as a single peak line centered at 5.7kms“^, with a measured 
FWHM of 1.9kms“^, and a maximum of 72.8mJy, based on a Gaussian fit to the line 
profile. The emission peaks from both the continuum and the integrated spectral line appear 
separated by 0.38 ± O.lOarcsec (Fig. [^. Integrating our line profile over the corresponding 
channels we obtain a total flux of 156mJykms“^. 


3. Spectral Energy Distributions 


The detection of gas around the third object in FWTau is consistent with the proposed 
accreting nature of the disk and provides additional information to potentially constrain the 
properties of the third object. Before investigating this, we here summarize the observational 
data available for the triple system FWTau with special emphasis on separating resolved 
and unresolved observations. 

The SED of the full FW Tan system is dominated by emission from the photospheres of 
the central binary at wavelengths shorter than ~15 /xm. However, high resolution observa¬ 
tions have also been able to identify the contribution from the third object in the J, H, K, 


and L bands (Kraus et al. 2014). Emission at 1.3 mm is also confirmed as originating from 
the third component, both in this paper and in previous observations (K2015). This study 


also And a lack of 1.3 mm flux from the central binary, allowing upper-limits to be derived 
for this component. 

The excess emission at wavelengths longer than 15/xm and shorter than 1.3 mm is un¬ 
resolved, i.e. it could stem from a circumbinary disk around the central binary and/or a 
disk around the third object. These unresolved observations cover the Spitzer photometry 


from Cieza et al. (2012), the Herschel data from Howard et al. (2013), and the 450 /xm upper 


limit and 850 /xm detection from Andrews & Williams (2005). Given the 1.3mm detections 


of the third object, it is very likely that the later detection is also associated with it. Table 
[T] summarizes the photometric data available for the triple system FW Tan separating their 
different contributions. 
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Disk models 


The first attempts to understand the nature of the third object in FW Tan object rapidly 


led to the hypothesis of it being a PMC (Kraus et ah 2014). However, as discussed in detail 


in 


K2015| it might also be a late-type stellar or sub-stellar companion embedded in an edge- 


on disk. We here confront both scenarios with the complete set of available observations 
by presenting an example model for each case. A complete parameter study is beyond the 
scope of this letter and not warranted given the limited quality of the current data. 

For the SED analysis we assumed an extinction towards the line of sight of Ay = 


0.4 (Kraus et al. 2014). We applied this correction to both the binary and third object 


photometry. The binary contribution was assumed to be the sum of two identical M5 dwarfs 
normalized to the J-band flux, with temperatures of Teg = 3200 K, masses of M = 0.22 Mq 
and a distance of 140pc. These assumptions imply both stars have a radius of 1.1 R©, 
suggesting an age of ~ 2Myr (Baraffe et al. 2015), which is in agreement with the Taurus 
age. 


4.1. An edge-on disk around a substellar object 


Assuming all the excess emission comes from the third object, the SED resembles that of 
an edge on disk around a (sub)stellar body. To model this possibility, we removed from the 
observed excesses the contributions from the central binary emission. Using the radiative 
transfer code MCFOST (Pinte et al. 2006, 2009) we hnd that for inclinations i > 85 deg 
the extinction on the central object is too high to agree with a spectral type later than M5 
as required from spectroscopy (Bowler et al. 2014). Therefore, we focus on models with 
inclinations in the range of 70-85 deg, in which the disk has a direct effect on the central 
object emission as seen from our line of sight without causing too much extinction. 


We explore various parameters for reproducing the SED, i.e. the disk dust mass Md, 
the scale height at 100 au Hq, haring angle /3, grain size distribution (amin, Omax) for silicate 
composition, and the surface density exponent 7 for a power-law disk geometry. We also 
explore diherent values for the temperature and luminosity of the central object. We hxed 
the size of the disk to be 100 au in diameter, in agreement with the maximum allowed by 
the non-resolved ALMA detections. The continuum emission is only weakly dependent on 
this assumption as the dust is optically thick only to a few tens of au, thus a smaller disk 
will still remain as a plausible alternative with small ehects on the model SED. However, 
the gas emission is optically thick throughout, and will be ahected by the assumed disk size. 
For all high inclination conhgurations we identify a degeneracy between the scale height. 
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flaring exponent, inclination, and central object mass, which makes the modeling of this disk 
particularly challenging. Nonetheless, we hnd that for most parameter combinations a dust 
mass of ~ 3 Mq and a disk inner radius of ~ 1 au are required. Similarly, we hnd that an 
inclination close to i ~ 80.5 deg is required, as slightly lower inclinations expose too much 
stellar hux and a slightly higher values extinct too much, making the model inconsistent 
with the observed SED. Although the mass of the central object is less constrained by the 
SED, it is possible to hnd reasonable representations of the SED for masses < 0.14 Mq. 


However, we can constrain the mass of the central object in the edge-on model if we 
require its age to be close to the one estimated for the central binary. The larger the assumed 
mass for the central object, the older it must be to reach agreement with the observations. 
Figure]^ shows this dependence, assuming the evolutionary models of Barahe et ah (2015) 
and the calculated photospheric luminosity of 9 x 10“^ Lq. Even a conservative upper-limit 
of 6 Myr (~ 3 times the estimated age of the central binary) would imply the mass of the 
third object to be < 60 Mjup. We can therefore rule out stellar masses but the third object 
might be a forming brown dwarf seen edge-on. 


Baraffe et ah (2015 


The detected (2-1) emission might provide a further observational constraint on 
the nature of the third body. Therefore, we created an MCFOST radiative transfer model 
for the detected gas emission, assuming an ISM CO/H 2 abundance, a freezing temperature 
of 20 K (Qi et al. 2004; de Gregorio-Monsalvo et al. 2013[ ), turbulent velocities ranging from 
0.1-0.2kms“^, and the same spectral resolution as provided by the observations. We found 
that, given the low signal to noise of the line and the effects of photodissociation by the 
third object and the central binary ( Williams fc Best||2014 ), the total gas content of the disk 
remains an observationally unconstrained parameter. 


Moreover, the shape of the line does not provide unambiguous observational constraints 
on the mass of the central object. For the spectral resolution given by our observations, 
masses > 35 Mjup predict double-peaked lines in disagreement with the observed single- 
peaked line. However, we can not exclude such masses as the high noise level of the detection 
could have smeared out the line splitting. In addition, the observed line emission could be 
signihcantly affected by cloud contamination. Indeed, given that the measured velocity of 
the third object in FWTau is close to the mean velocity of the Taurus cloud, i.e. ~ 7kms“^ 
(Goldsmith et al. 2008), it is possible that either the blueshifted or redshifted emission from 
the disk and the cloud appear fused in the interferometer, thus precluding the interferometer 


from resolving the individual contributions and hltering them out (see e.g. Ganovas et al. 


2015). The same effect could explain the small (2.4 (t) difference in position of the centers of 


the continuum and line detections (see 


As an example for an edge-on model that explains all the available observations, we 






















7 


show a model fit to the SED and the CO line assnming a 5Myr old 40Mjup third object 
in Fig.|^ where it is clear that all the observed excess flnx can be explained by material 
confined to the third object. 


4.2. A disk around a planetary mass companion object 


The third object in FWTan conld also be a cooler low-mass object with a disk seen at 
low inclination. Indeed, the optical and near-IR emission is well fitted by a Tgff ~ 1800 K 
snbstellar model, with a Inminosity implying a 1-2 Myr old age, in agreement with the age 


of Tanrns forming objects, and consistent with a 5-7Mjup mass body (Baraffe et al. 2015). 


The projected separation of the third object from the binary has been measnred to be 
330 an, and if the trne separation is close to the projected one, the flnx from the nearby binary 
will contribnte to disk heating, thns affecting the temperatnre strnctnre and conseqnently 
its emission. To reprodnce the observed SED we have modeled the disk aronnd the third 
component assnming the binary contribntion to be identical to that described in 14.1 The 


separation between the binary stars was assnmed to be 11 an, whereas the distance and the 
relative orientation of the binary with respect to the disk were changed in order to explore 
different conhgnrations for the external illnmination. Apart from this additional parameter 
we explored the same parameter space as in 


1.1 


We hnd that the SED of FWTan can be explained with models in which the binary 
illnminates the disk with an incident angle ~ 20 — 50 deg with respect to the disk plane as 
these orientations maximize the heating of the disk (Fig. [^. Other conhgnrations do not 
expose the disk to enongh radiation, creating a lack of emission at longer wavelengths. We 
hnd that ~ 2 — 3 M 0 of dnst are enongh to reprodnce the emission at (snb) mm-wavelengths. 


in agreement with the estimates of K2015 The mid/far-IR emission is consistent with very 
small haring angle exponents (i.e. /3 ~ 1.0), and an inner radins of ~ 0.3 an. The overall hnx 
excess is only achieved, however, by introdncing a relatively large scale height (~ 22 — 25 an 
at 100an), and models with smaller valnes tend to nnderestimate the mid/far-IR emission. 

Concerning the gas line, for an assnmed central body mass of 7Mjup, the spectral 
prohle appears as a clear single line for inclinations i < 15 deg, in agreement with the 
observations (see Fig.[^, and the Keplerian donble-peak appears only for larger inclinations. 
The amplitnde of the CO line peak snggests a very small gas mass, becanse the amplitnde of 
the line predicted by the model, for an ISM CO/II 2 abnndance and gas-to-dnst mass ratio, is 
mnch larger than the observed valnes. This conld be either an ehect of CO photodissociation 


which can rednce the CO/II 2 abnndance (Visser et al. 2009 Williams & Best 2014), an 

















intrinsic very low gas-to-dust mass ratio, and/or a smaller disk size because a smaller emitting 
area would predict less CO flux and thus require a less reduced gas content. Similar to the 
high-inclination case, it is also possible that the lack of flux from the observations compared 
to the model might be due to the interferometer Altering out a major fraction of the flux. 


5. Conclusion 

We present ALMA Cycle 1 band 6 continuum and ^^CO (2-1) observations of the triple 
system FWTau and test two different scenarios against the new observations. 

We And that the edge-on model can recreate all observed excess emission. Although 
the observed single-peaked line is most naturally reproduced assuming a low mass for the 
third object, larger masses can not be excluded because of the low SNR of the detection and 
the possibility of cloud contamination. Assuming the third object and its central binary are 
coeval however, evolutionary tracks strongly suggest it is a substellar object. 

The low inclination scenario for a planetary mass object, in combination with the radia¬ 
tion received from the close binary, is capable of reproducing the SED and the single-peaked 
line emission. The weakness of the (2-1) line either indicates a low gas mass, significant 
photodissociation affecting the CO/H 2 abundance, significant contamination by the cloud, 
a smaller disk, or a combination of these four effects are at work. 

Resolved continuum observations of both the binary and third components of the 
FW Tan system, as well as resolved optically thin line detections including both CO higher- 
transition and isotopologues, will definitively unveil the nature of this intriguing system. 
Regardless, all data suggest FWTau is a substellar object caught in formation, cementing 
its position as a vital object for understanding the brown dwarf/planet formation process. 
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erated under cooperative agreement by Associated Universities Inc. This publication makes 
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Table 1. FWTau photometry data 


Wavelength 

(/im) 

Flux 

FW Tan A/B 

Flux 

Third component 

Flux 

Unresolved 

Unit 

Reference® 

1.25 

10.34 

17.34 ±0.07 


mag 

1 

1.66 

9.67 

16.24 ±0.07 


mag 

1 

2.20 

9.38 

15.29 ±0.07 


mag 

1 

3.35 



9.196 ±0.026 

mag 

2 

3.6 



64.8 ±0.3 

mJy 

3 

3.8 

9.19 

14.25 ±0.10 


mag 

1 

4.5 



44.8 ±0.1 

mJy 

3 

4.6 



41.5 ±1.0 

mJy 

2 

5.8 



32.4 ±0.2 

mJy 

3 

8.0 



18.0±0.1 

mJy 

3 

12 



8.2±0.3 

mJy 

2 

22 



< 10.5 

mJy 

2 

24 



6.79 ±0.4 

mJy 

4 

70 



30 ±4 

mJy 

5 

100 



33 ±4 

mJy 

5 

160 



70 ±40 

mJy 

5 

450 



< 30.0 

mJy 

6 

850 



4.5 ±1.1 

mJy 

6 

1300 

< 0.084 

1.78 ±0.03 


mJy 

7 

1300 


3.4±0.2 


mJy 

This work 


Note. — 

“References 

(1) Kraus et ah ( 

2014); (2) Wright et ah 

(2010); (3) Cieza 

et ah (2009); (4) Cieza et al.|(2012); (5) Howard et al.](2013); (6) Andrews & Williams 

(2005|); (7) Kraus et ah 

(2015 

)• 
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Table 2. Parameters in the models 


Parameter 

PMC model 

Edge-on disk model 

Unit 

i 

10 

80.5 

deg 

Ms 

7 

40 

Mjup 

dPdisk 

3 

3 

M© 

Ho 

25 

15 

an 


1.00 

1.25 


7 

-0.65 

-0.1 


^min 

1.0 X 10-2 

1.0 X 10-2 

/im 

^max 

1.0 X 10^ 

1.0 X 10^ 

/im 


Note. — These are the parameters used to create the 
models shown in Fig. 



3 2 1 0 -1 -2 -3 

Relative right ascension (arcsec) 


Fig. 1.— Integrated ^^CO(2-l) emission map. The solid-line contours show the continuum 
detection above 5a in steps of 2a. The dashed contour limits the 3a emission from the (2-1) 
line used to calculate the line profile. The white ellipse shows the synthesized ALMA beam. 
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Fig. 2.— Age-mass region constrained by the edge-on modeling. Gray regions represent the 
1 — 3cj range of allowed ages for a given central object mass. The red-dashed line indicates the 
deuterium-burning limit mass for a solar composition (Chabrier et al. 2000). The shaded area 
limits the maximum allowed mass given spectroscopic constraints (Bowler et al.|pM4). 
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Fig. 3.— Left panels: SED for the FWTau system. Green squares represent the emission from the 
central binary, blue squares represent the emission associated to the third component, cyan circles 
and triangle are the unresolved photometric observations and upper limit respectively. Upper-left 
panel shows the edge-on model described in 14.1 (blue solid line), including the highly absorbed 
substellar photosphere (gray solid line). Lower-left panel shows the PMC model (red solid line) 
including the photospheric emission from a 1800 K central object, and the integrated SED model 
composed of the PMC model plus the photospheric emission from the central binary, shown as 
M5-type photospheres (^.2). The parameters used to create these models are specified in Table 
Right panels: Detected ALMA ^^CO (2-1) line profile measured over the 3a region shown in Fig. [b 
(black solid line). The dotted lines represent the edge-on model ( ^4.1[ upper-right panel), and the 
planetary mass companion model ( ^4.2[ lower-right panel), respectively. The unconstrained total 
gas mass was adjusted to fit the maximum detected emission. 
































